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Long-term hematopoietic stem cells (HSCs [LT-HSCs]) are well known to display unpredict- 
able differences in their clonal expansion capacities after transplantation. Here, by analyz- 
ing the cellular output after transplantation of stem cells differing in surface expression 
levels of the Kit receptor, we show that LT-HSCs can be systematically subdivided into two 
subtypes with distinct reconstitution behavior. LT-HSCs expressing intermediate levels of 
Kit receptor (Kit int ) are quiescent in situ but proliferate extensively after transplantation 
and therefore repopulate large parts of the recipient's hematopoietic system. In contrast, 
metabolically active Kit hi LT-HSCs display more limited expansion capacities and show 
reduced but robust levels of repopulation after transfer. Transplantation into secondary and 
tertiary recipient mice show maintenance of efficient repopulation capacities of Kit int but 
not of Kit hi LT-HSCs. Initiation of differentiation is marked by the transit from Kit int to Kit hi 
HSCs, both of which precede any other known stem cell population. 



Hematopoietic stem cells (HSCs) replenish 
millions of mature hematopoietic cell types 
every second throughout life but also maintain 
the HSC pool over time. HSC function is as- 
sessed by their capacity to repopulate the blood 
system of lethally irradiated recipient mice in 
the long term. The most immature HSC pool is 
functionally heterogeneous, and HSCs vary in 
their differentiation potential and duration of 
reconstitution (Copley et al., 2012; Muller- 
Sieburg et al., 2012). However, the magnitude 
of repopulation, thus white blood cell output 
per donor HSC, was only retrospectively asso- 
ciated with specific reconstitution patterns de- 
termined by lineage choice (Dykstra et al. , 2007) . 
Therefore, it remains unknown whether clonal 
expansion capacities are predetermined in donor 
cells or whether the magnitude of repopula- 
tion is determined by the microeiivironment of 
the recipient. 

Kit expression is widely used for the pro- 
spective isolation of HSCs, and the stem cell 
factor (SCF)— Kit signaling axis is pivotal for 
normal pool size and function of fetal and adult 
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HSCs (Russell, 1979; Ikuta andWeissman, 1992). 
Consistently, alterations in Kit signaling pro- 
foundly affect adult HSC function (Ogawa 
et al., 1991; Czechowicz et al., 2007; Waskow 
et al, 2009; Ding et al., 2012; Deshpande et al, 
2013). Furthermore, Kit alleles resulting in hy- 
pomorphic expression of the receptor are loss 
of function alleles (Russell, 1979;Thoren et al., 
2008; Waskow et al., 2009), suggesting that re- 
duced densities of Kit expression correlate with 
loss of "sternness." In contrast, cells expressing 
low levels of (Doi et al., 1997; Matsuoka et al., 
2011) or lacking (Ortiz et al., 1999) Kit recep- 
tor expression were suggested to contain quies- 
cent long-term HSCs (LT-HSCs). However, 
differences in the clonal expansion capacities of 
HSCs expressing distinct levels of the Kit re- 
ceptor were not reported. 
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Figure 1. Prospective separation of LT-HSCs with high and low clonal expansion capacities based on the level of Kit expression. (A) BM was 

harvested from B6 mice, and Kit int and Kit hi cells within the LSK Slam population were sorted based on the indicated gates. (B-D) 50 sorted LSK Slam Kit int 
or Kit hl cells were transplanted into primary (1°), secondary (2°), and tertiary (3°) recipient mice, and donor cell contribution to the blood neutrophil com- 
partment (B), HSC compartment (C), and donor cell composition in the blood (D) were analyzed by flow cytometry. Number of recipient mice per donor 
population: 1°, 3; 2°, 6; and 3°, 6. T cells (CD3+), B cells (B220+), and myeloid cells (Myel.; CD 1 1 b + ) are shown in D. n = 3 per condition. Mean + SD is rep- 
resentative of two independent experiments. (E) BM cells from donor mice (WT Ctrl), post sort cells before transplantation (open squares), and HSC progeny 
in recipient mice (closed circles) were isolated and analyzed for the expression density of Kit (MFI). Expt 1 (top): 50 LSK CD48~ CD41 " CD150 + Kit"" or Kit 1 " 
donor cells. Expt 2 (bottom): 1,000 LSK Kit int or Kit hi donor cells. Error bars show SD. (F) LSK Kit int or LSK Kit hi cells were sorted from congenic wild-type 
mice, mixed, and transplanted into irradiated recipients. 18 h later, contribution of Kit int - or Kit hi -derived cells within all donor LSK cells in the recipient's 
BM was determined. Mean + SD of three recipient mice is shown. Data are representative of two independent experiments. (G) LSK Slam Kit int or Kit hi cells 
were sorted, and limited cell numbers (top: 3, 8, and 20; bottom: 3, 8, and 13) of each population were transplanted together with 3 x 10 s congenic BM 
cells into 10 lethally irradiated recipients each. Contribution of test donor cells to neutrophils was analyzed 16 wk later. Data are from two independent 
experiments. *, P < 0.05; **, P < 0.01 ; ***, P < 0.001. 



To assess whether expansion capacities are predetermined 
within donor HSCs and whether this function identifies 
novel cellular subtypes within the most immature HSC pool, 
we transplanted LT-HSCs that differed in the density of the 
expression of the Kit receptor. Donor cells repopulated re- 
cipient mice to two significantly different magnitudes: HSCs 
with intermediate levels of Kit receptor expression (Kit mt ) 
contained greater expansion capacities compared with HSCs 
expressing high densities of the Kit receptor (Kit hl ), suggest- 
ing that HSC clonal growth potential is predetermined in a 
cell-intrinsic fashion. We further provide evidence that these 
HSC subtypes are two consecutive developmental stem cell 
stages within the most immature HSC pool and that transit 
from Kit mt to Kit 111 LT-HSCs marks the onset of differentia- 
tion and is associated with significant loss of expansion capac- 
ities. Gene expression profiles ex vivo and after SCF trigger 
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suggest that the inherent differences are based on distinct cy- 
cling and adhesive activities. 

RESULTS AND DISCUSSION 

Prospective separation of HSCs with different expansion 
capacities: Intermediate levels of Kit receptor expression 
correlate with increased HSC potency 

To assess whether distinct levels of Kit cell surface expression 
mark discrete types of HSCs that differ in their biological 
properties, we fractionated the HSC compartment into cells 
expressing high and intermediate densities of the Kit recep- 
tor (Fig. 1 A) and performed competitive transplantation 
experiments. Both donor populations engrafted stably over 
time (Fig. 1 B). However, Kit mt cells showed high repopula- 
tion of blood neutrophils and BM-resident HSCs, whereas 
Kit' 11 cells contributed to sustained but low levels in both 
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Figure 2. Density of Kit receptor expression defines two functionally distinct populations. (A) BM cells were prepared, and expression of the indi- 
cated markers was determined on Kit mt and Kit hl cells pregated as described above the plots. Isotype control staining on LSK Slam cells is shown in black. Data 
are representative of two independent experiments each. (B) BrdU was injected into wild-type mice, and after 4 h, BM cells were prepared and BrdU incorpo- 
ration into LS Kit + cells or into LS Kit lnt and Kit 1 " cells was measured. (C) Plot shows summary of four mice analyzed in two independent experiments as out- 
lined in B. (D) LSK cells were labeled with BrdU, and 330 d later, Kit expression on LS CD48~ CD150+ BrdU-positive and -negative cells was analyzed. Data are 
representative of two independent experiments using three mice each. (E) Lineage-depleted BM cells were cultivated with or without SCF for 1 5 min, and 
phosphorylation of Akt (top) and Erk (bottom) in LSK Slam Kit" 11 (left) or Kit hl (right) HSCs was measured. (F) Plot summarizes data from three independent 
experiments using cells from 10 (pErk) or 1 1 (pAkt) mice as shown in E. For controls, four (pErk) or five (pAkt) mice were used. (G) LSK CD135~ Kit lnt or Kit hl 
HSCs were cultivated overnight in the presence or absence of SCF, BrdU was added for 4 h, and its incorporation was subsequently analyzed. Specificity was 
determined by using the anti-Kit inhibitory antibody ACK2 (Ogawa et al., 1991). The plot summarizes data from two independent experiments using three 
mice for each group. (H) Single LSK Slam CD135~ CD34~ Kit + , Kit"", or Kit hl HSCs were cultivated in liquid culture in medium supplemented with SCF, Tpo, 113, 
and Epo (left). Colony formation was analyzed after 14 d. L, large; M, medium; S, small colonies. Cell type composition of individual colonies is summarized 
(right), m, macrophage; g, granulocyte; E, erythroblast; M, megakaryocyte. Data shown are pooled from two independent experiments using cells from three 
mice each. (I) LSK CD135~ Kit lnt or Kit hl HSCs were cultivated as described in G, and the frequencies of Kit + Sca-1 + -expressing cells were determined. Plot 
summarizes data from three mice. Changes are quantified in the plot on the right. Error bars show SD. *, P < 0.05; **, P < 0.01 ; ***, P < 0.001. 



compartments (Fig. 1 C). Donor cell contribution was stable 
for Kit 11 " HSCs and their progeny in secondary and tertiary 
recipients, whereas contributions of Kit hl -derived HSCs 
declined over time and eventually became nondetectable. There 
was no difference in the composition of Kit 1111 - or Kit 1 "- 
derived mature white blood cells (Fig. 1 D). Expression den- 
sities of the Kit receptor on donor HSCs derived from 
Kit 11 " HSCs (Kit lnt -Sort, Fig. 1 E) were increased in pri- 
mary recipient mice 16 wk after transplantation and subse- 
quently gradually declined on HSC progeny in secondary 
and tertiary recipient mice. In contrast, Kit expression levels 
on donor HSCs derived from Kit hl HSCs (Kit hl -Sort) were 



already declined in primary recipient mice, suggesting that 
Kit 111 ' HSCs precede Kit 1 " HSCs during differentiation. 

Disparate repopulation activities after bulk transplantation 
can be caused by different homing activities or by varying 
HSC frequencies within the transplanted populations. How- 
ever, homing activity between Kit"" and Kit 1 " cells was undis- 
tinguishable (Fig. 1 F), and HSC frequency was also comparable 
as determined by limited dilution transplantations (Fig. 1 G; Hu 
and Smyth, 2009), suggesting that differences in repopulation 
activities are caused by cell-intrinsic differences in the clonal 
expansion capacities of donor cells. Collectively, these results 
indicate that clonal expansion capacities are cell-intrinsically 
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Figure 3. Kit receptor densities are indicative for HSC populations with different molecular expression programs. BM was harvested and Kit "' 
or Kit hi LT-HSCs (LSK CD48~ CD41" CD150+ CD135- CD34"), ST-HSCs (LSK CD48+ CD135- CD34+), and MPP cells (LSK CD48+ CD135+ CD34+) were sorted. 
Whole transcriptome sequencing was performed on amplified RNA; n = 3 for each population. (A) Heat map of the sample to sample Euclidean distance 
and dendrogram showing sample to sample correlation. All biological replicates cluster well with each other, and all samples of the different populations 
are clearly separated from each other. (B) Graph shows DESeq-normalized counts for selected genes expressed in Kit lnt or Kit hl LT-HSCs. Protein expression 
from the same genes is depicted in Fig. 2 A. Mean ± SD is shown; n = 3; **, P < 0.01. (C) Comparison of gene expression profiles between Kit' nt and Kit hi 
LT-HSCs. Red dots represent significant DEGs at a 10% false discovery rate (FDR). (D and E) Plots show biological processes that are enriched in genes 
down-regulated (D) or up-regulated (E) in LT-HSC Kit lnt cells compared with LT-HSC Kit hi cells. Analysis was performed using the GO/BP database of DAVID. 
Enrichment scores (—log transformation of the DAVID EASE score) were calculated to determine overrepresentation of particular biological processes 
and are indicated on the x axis. (F) Plot shows biological processes that are uniquely enriched in genes differentially expressed in LT-HSC Kit int cells after 
stimulation with SCF compared with LT-HSC Kit 1 " cells after stimulation with SCF. (G) Heat map shows quantitative comparison of the expression level of 
genes common between DEGs of Kit lnt ± SCF and Kit hl ± SCF that are associated with the top five enriched GO biological process terms for these common 
genes. The number of associated DEGs is indicated in parentheses. 



predetermined and that different expression densities of the 
Kit receptor allow for the prospective separation of two sub- 
types of LT-HSCs differing in that function. 

Level of Kit expression marks two 
functionally distinct HSC populations 

HSC subtypes with prolonged self-renewal activity can be 
found in stem cell populations that express low levels of 
CD49b (Benveniste et al, 2010) or high levels of CD150 
(Morita et al, 2010), CD86 (Shimazu et al, 2012), and 
CD41 (Gekas and Graf, 2013). Furthermore, HSCs that are 
biased toward the development of either lymphoid or my- 
eloid cells respond differently to stimulation with IFN-ct 
(Essers et al., 2009) or TGF-(3 (Challen et al, 2010). The 
latter can modulate Kit cell surface expression (Sansilvestri 
et al., 1995). However, expression of receptors for TGF-(3 
(Tgfbr2) and IFN-a (Ifnarl) was comparable on Kit mt and 
Kit' 11 HSCs (Fig. 2 A). CD 150 expression was inversely cor- 
related with Kit expression, but there was no evidence for 



differential expression of CD49b, CD41, CD86, or other 
stem cell markers like Epcr (Proa, CD201) and Seal (Ly6a) 
on LSK Slam Kit mt or Kit hl cells, suggesting lack of overlap 
between these populations. 

Next, we showed that Kit"" cells were largely quiescent 
and in the G0/G1 phase of the cell cycle and that only Kit hl 
cells contained actively cycling cells (Fig. 2, B and C). Con- 
sistent with this cell cycle profile, BrdU label— retaining cells 
were found to express low densities of the Kit receptor after 
a chase period of 330 d (Fig. 2 D), suggesting that Kit mt cells 
rarely cycle under steady-state conditions. We conclude that 
differential expression of Kit allows for the separation of ac- 
tive and quiescent HSCs. 

Intensity of signal transduction depends 
on the density of Kit receptor expression 

To analyze whether distinct cell biological properties may be 
caused by differential SCF— Kit signaling between Kit lnt and Kit' 11 
cells, we determined the frequency and the mean fluorescence 
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intensity (MFI) of phosphorylated Erk and Akt after SCF 
trigger (Fig. 2 E).The frequency of Kit"" cells that phosphory- 
lated Erk after SCF triggering was reduced compared with 
Kit hl cells, and furthermore, the MFI of pErk and pAkt was 
decreased in Kit lnt cells, suggesting that reduced densities of 
Kit cell surface expression result in decreased signaling activ- 
ity (Fig. 2 F). Also, cell biological consequences after SCF— 
Kit signaling were found to be different between both 
populations: Kit"" cells proliferated to lower rates after stimu- 
lation with SCF compared with Kit'" cells (Fig. 2 G), and a 
reduced frequency of Kit"" cells gave rise to colonies after 
14 d of culture (Fig. 2 H). Differential signaling and proliferation 
depended on SCF-mediated effects because both activities 
were abrogated by the inhibition of Kit signaling using a 
pharmacological inhibitor (Fig. 2 F) or a blocking anti-Kit 
antibody (Fig. 2 G). Finally, reduced signaling activity in 
Kit 11 " cells inversely correlated with the maintenance of a Kit + 
Sca-1 + stem cell surface phenotype after culture (Fig. 2 I), 
suggesting that high levels of Kit signaling result in differenti- 
ation but low levels of Kit signaling result in maintenance of 
sternness. We conclude that different densities of Kit receptor 
expression result in altered cell biological consequences after 
triggering the SCF— Kit signaling axis. 

Kit int and Kit hi LT-HSCs are distinct molecular entities 

To test for molecular differences between both HSC popula- 
tions, we compared the gene expression between Kit"" and Kit 1 " 
LT-HSCs, short-term HSCs (ST-HSCs), and multipotent pro- 
genitor (MPP) cells (Fig. 3). Unsupervised clustering revealed 
great homology between individual samples (Fig. 3 A), and 
analysis of expression counts of selected transcripts encoding for 
defined proteins that were previously analyzed by flow cytom- 
etry (Fig. 2 C) confirmed the purification strategy of both pop- 
ulations (Fig. 3 B). Comparing gene expression between Kit"" 
and Kit'" HSCs revealed 96 up-regulated and 48 down-regulated 
genes in Kit 1 "' HSCs (Fig. 3 C).The enriched functional an- 
notation terms associated with genes differentially expressed 
between both populations showed an overrepresentation of 
terms that map to pathways related to cell cycle/division for the 
significandy down-regulated genes, verifying the quiescent state 
of HSCs that exhibit great expansion capacity after transplanta- 
tion (Fig. 3D). In contrast, we found an overrepresentation of 
terms related to cell adhesion pathways for the significandy up- 
regulated genes (Fig. 3 E), supporting the idea that retention of 
HSCs in their niche space instructs quiescence. 

To assess whether signaling via the SCF— Kit axis results in 
different gene expression responses in Kit"" or Kit 1 " HSCs, we 
compared differentially expressed genes (DEGs) between both 
cell types after stimulation with SCF. Consistent with the in- 
duction of cell division upon SCF trigger, DEGs that were 
uniquely found in Kit"" HSCs showed an overrepresentation 
of terms that map to pathways related to cell cycle (Fig. 3 F). 
Furthermore, DEGs common to Kit"" and Kit 1 " HSCs showed 
a reduced expression score in Kit"" cells, indicating that the 
genes that are induced in both populations are expressed at 
lower levels in Kit"" cells (Fig. 3 G). This finding is consistent 
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Figure 4. Present and modified model for the initiation of HSC 
differentiation based on transplantation experiments. (A) Present 
model: Kit-positive (plus) LT-HSCs repopulate recipient mice in the long 
term (each tick on the x axis represents ~1 6 wk). Kit+ intermediate-term 
reconstituting cells (ITRCs), ST-HSCs, and MPP cells repopulate recipient 
mice to lower levels and for much shorter time periods than LT-HSCs. The 
repopulation scheme of LT-HSCs, ST-HSCs, and MPP cells is based on the 
transplantation of 1 ,500 cells per population (not depicted). The repopula- 
tion scheme of ITRCs was estimated from Benveniste et al. (2010). (B) Mod- 
ified model: We show that the magnitude and duration of reconstitution 
differs between HSC subsets and that the most immature LT-HSC compart- 
ment can be further subdivided into two novel stem cell subtypes that have 
inherent differences in their clonal expansion capacities. Kit"" HSCs re- 
populate recipient mice very efficiently (magnitude and duration), whereas 
Kit 1 " cells can be placed between Kit"" LT-HSCs and ITRCs. Both magnitude 
and duration of reconstitution are reduced compared with Kit'"' LT-HSCs 
but increased compared with ITRCs. Q, quiescent LT-HSC. 

with lower levels of Kit signaling in Kit"" HSCs directly after 
trigger with SCF compared with Kit 1 " HSCs and suggests that 
differential gene expression in both populations is a direct 
consequence of different levels of Kit signaling in situ. 

Collectively, we identified two discrete subtypes of HSCs 
that differ in their inherent capacities to expand and to form 
HSC progeny after transplantation. Kit"" HSCs are the most 
immature HSCs that may be parental to Kit' 11 HSCs. We 
show that Kit expression levels negatively correlate with 
sternness and that reduced levels of signaling via the SCF— Kit 
axis are required for the maintenance of HSC function. Tran- 
sition between both HSC subtypes is accompanied by a dis- 
tinct loss of repopulation potential based on different clonal 
expansion capacities (Fig. 4). Both HSC subtypes are parental 
to any previously described HSC population, suggesting the 
identification of a novel cell type at the top of the hemato- 
poietic hierarchy that marks initiation of differentiation. The 
separability of both populations will help to understand the 
mechanisms of self-renewal and it will be interesting to de- 
termine whether the frequency and function of both LT- 
HSC subsets remain constant over time during aging. 
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MATERIALS AND METHODS 

Mice. C57BL/6 (B6) and B6.SJL-Ptprt'Pep3 h /Boyj (B6.SJL) mice were pur- 
chased from the Jackson Laboratory and bred and maintained under specific 
pathogen— free conditions in the animal facility at the Medical Theoretical 
Center of the University of Technology Dresden. Experiments were per- 
formed in accordance with German animal welfare legislation and were ap- 
proved by the relevant authorities, the Landesdirektion Dresden. 

Transplantation. For competitive transplantation, 50 lineage" Sca-1 + Kit + 
(LSK) CD48" CD41" CD150+ (Slam) Kit?* or Kit hl cells (Fig. 1, B-E [top]) 
or 1,000 purified LSK Kit 11 " 1 or Kit' 11 cells (Fig. 1 E, bottom) were transplanted 
together with 5 x 10^ no nfrac donated BM cells into lethally irradiated (900 
cGy) wild-type recipients. Test, competitor, and recipient cells carried differ- 
ent CD45 alleles (CD45.1 + , CD45.2 + , CD45.1 + CD45.2+). For serial trans- 
plantation, 5 x 10 6 nonseparated BM cells were injected into secondary 
lethally irradiated recipients. For limiting dilution analysis, 3, 8, and 20 or 3, 
8, and 13 LSK Slam Kit im or Kit" (CD45.1 + ) cells were injected together 
with 3 x 10 5 nonseparated BM cells (CD45.1 + CD45.2 + ) into 10 lethally ir- 
radiated wild-type mice (CD45.2 + ) per donor cell number. 16 wk after trans- 
plantation, donor cell chimerism was determined in blood neutrophils, and 
mice were scored positive when donor contribution was >1%. Frequency of 
the repopulating cells was calculated using ELDA software. Pairwise differ- 
ences in active cell frequencies between groups were calculated as described 
previously (Hu and Smyth, 2009). For the first experiment, donor cells were 
sorted into one well and separated before transplantation, and for the second 
experiment, donor cells were sorted into separate wells of a 96-well plate. 
For homing assays, LSK Kit™ 1 (CD45.1 + or CD45.2+) and LSK Kit hi 
(CD45.2* or CD45.1 + ) cells were sorted and mixed. Ratio of mixture was 
determined by flow cytometry, and a total of 2 X 10 4 cells were injected into 
each lethally irradiated recipient mouse (CD45.1 + CD45.2 + ). 16—18 h later, 
recipients were sacrificed, BM cell suspension was prepared, and donor cell 
ratio in LSK cells was determined. 

Flow cytometry. BM cell suspensions were prepared, stained, and analyzed 
as described previously (Arndt et al., 2013). Antibodies (clones in parenthe- 
ses) used are as follows: CD3 (2C11; 17A2), CD lib (Ml/70), CD19 (1D3), 
CD34 (RAM34), CD45.1 (A20), CD45.2 (104), CD45R (RA3-6B2), 
CD86 (GL1), CD117 (2B8), CD135 (A2F10), Gr-1 (RB6-8C5), Nkl.l 
(PK136), Seal (D7), Terll9 (Terll9), CD41 (MWReg30), Epcr (1550), 
CD49b <DX5), and Ifnrl (MAR1-5A3; all eBioscience); CD48 (HM48-1) 
and CD150 (TC15-12F1; BioLegend); and Tgfbr2 (polyclonal; R&D Sys- 
tems). MFI of Kit receptor expression was normalized between independent 
experiments based on the MFI of Kit expression on wild-type LS CD48" 
CD150 + cells: {MFI Kit on donor-derived cells in experiment 1) = (MFI Kit 
on wild-type cells in experiment 1)/(MFI Kit on wild-type cells in experi- 
ment 2) x {MFI Kit on donor-derived cells in experiment 2). 

BrdU labeling in vivo. Mice were intraperitoneally injected with a single 
dose of BrdU (Sigma -Aldrich; 1 mg in 200 ul PBS) and sacrificed 4 h later. 
Cell cycle and BrdU incorporation analyses were performed as described 
previously (Waskow et al., 2008). To test for label retention, mice were in- 
jected once with BrdU as described and fed BrdU-containing drinking water 
(1 mg/ml). BrdU-containing water was replaced every third day, and after 
13 d replaced by normal drinking water. Mice were sacrificed after 330 d and 
BM cells were stained. A control group of three mice was sacrificed after 13 d, 
and 98 ± 2% of LSK cells had incorporated BrdU. 

In vitro culture. To analyze phosphorylation of Akt and Erk (BD), 2—4 X 10 6 
lineage-depleted BM cells/ml were cultivated for 15 niin in DMEM 
supplemented with 2% FCS and 50 uM [3-mercaptoethanol with or without 
100 ng/ml rmSCF (R&D Systems) and with or without 100 mM Imatinib 
(LC Laboratories) and subsequently analyzed as described by the manufac- 
turer. In brief, cells were fixed in prewarmed Lyse/Fix buffer (BD), washed 
and permeabilized using Perm Buffer III (BD), and stained for phosphorylated 
signaling molecules and cell surface markers for 20 min at room temperature. 



After a washing step, LSK Slam Kit mt or Kit' 11 cells were immediately analyzed 
by flow cytometry. To analyze BrdU incorporation, sorted LSK CD 135" 
Kit mt or Kit 1 " cells were cultivated overnight in StemSpan medium (STEM- 
CELL Technologies) supplemented with 100 ng/ml rmTpo, 100 ng/ml rmFlt3 
ligand, 50 ng/ml rmSCF, and controls with an additional 10 u.g/ml of an in- 
hibitory anti-Kit antibody (ACK2; eBioscience). Subsequently, 10 uM BrdU 
was added to the cultures for 4 h, and BrdU incorporation and cell surface 
phenotype were analyzed. For colony growth, single LSK Slam CD34" 
CD135" Kit im , Kit lu , or Kit + cells were sorted into individual wells of a 96-well 
plate into StemSpan medium supplemented with 20 ng/ml rmSCF, 20 ng/ml 
rmTpo, 20 ng/ml rmI13, and 5 U/ml rhEpo and cultivated for 14 d. Colony 
size was determined as follows: small (S) < 1 mm, medium (M) = 0.5—1 mm, 
and large (L) > 1.5 mm. Subsequently, cells from each colony were cyto- 
spun, and cell types were determined after May-Grunwald Giemsa staining as 
described previously (Arndt et al., 2013). 

RNA isolation, amplification, and sequencing. For RNA isolation, 
LSK Slam CD135" CD34" (LT-HSC) Kit™* or Kit hl cells, LSK CD48+ 
CD135" CD34 + (ST-HSC), and LSK CD48 + CD135 + CD34* (MPP) cells 
were sorted (FACSAria II; BD) and immediately lysed in jjlMACS mRNA 
isolation lysis buffer (Miltenyi Biotec). For SCF trigger, 7,500 sort-purified 
LSK CD135 - Kit im or Kit' 11 cells were incubated in StemSpan medium 
(STEMCELL Technologies) supplemented with 50 ng/ml rmSCF for 11 h 
and subsequently lysed. Lysates were cleaned using LysateClear Columns 
(Miltenyi Biotec), and mRNA was directly isolated from the lysis buffer 
using SeraMag oligo(dT14) beads (Thermo Fisher Scientific). The mRNA 
was eluted in a volume of 5 ul of 10 mMTris-HCL and directly subjected to 
subexponential RNA amplification using the WT-Ovation System (Nugen 
Technologies). Samples were prepared according to the manufacturer's in- 
structions, but stopped before final Post-SPIA Modification. After bead-based 
purification (XP beads; Agencourt), randomly primed second strand synthe- 
sis was performed using second Strand Synthesis Module from New England 
Biolabs, Inc. After DNA shearing by ultrasonication (Covaris S2) and treat- 
ment with SI nuclease (New England Biolabs, Inc.), samples were subjected 
to standard Dlumina fragment library preparation using indexed adaptors. 
Resulting libraries were pooled in equimolar quantities for 75-bp single-read 
sequencing on IUumina HiSeq 2000 and distributed on several lanes, result- 
ing in ^30—90 million reads per sample. 

Bio info rmatic analysis. Alignment of the short reads to the mm9 tran- 
scriptome was performed with pBWA software, and a table of read counts 
per gene was created based on the overlap of the uniquely mapped reads 
with the Ensembl Genes annotation version 61 for mm9, using BEDtools 
(version 2.11; Quinlan and Hall, 2010). The raw read counts were then 
normalized with the DESeq R package (version 1.8.1; Anders and Huber, 
2010), and the sample to sample Euclidean distance was computed based 
on the normalized counts to explore sample to sample correlation. After 
normalization, testing for differential expression was performed with 
DESeq, and accepting a maximum of 10% false discoveries (10% FDR) re- 
sulted in 96 up-regulated and 48 down-regulated genes in Kit im HSCs. For 
SCF trigger experiments, DEGs between Kit mt versus Kit mt + SCF and 
Kit 1 " versus Kit 1 ' 1 + SCF were compared, and DEGs unique to Kit ,m + SCF 
versus Kit mt + SCF or DEGs common to both gene lists were identified. 
To identify enrichment for particular biological processes associated with 
the DEGs, the DAVID GO/BP/FAT database (Huang et al., 2009) was 
used. Enrichment scores were calculated (— log transformation of the DAVID 
EASE score) to determine overrepresentation of particular biological pro- 
cesses. To quantify gene expression levels, an expression score defined as 
the median of all normalized counts of the DEGs associated with that par- 
ticular GO term was calculated. Subsequently, expression scores for the 
top five GO biological process terms across the four experimental condi- 
tions were depicted in a heat map. 
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